CARB  Mol  070*41  $t 


wwtti  1 2.  report  oate 

I  1991  March  22 


3.  WOIT  TWI  AMO  OATES  COVEAEO 

Annual  03/90  -  02/91 


a.  mu  and  suavmi 

Pressureless  Densificatictn  of  Ceramic  Harris  Composites 


B.  AUTHOMtS) 

M.  Ji.  Rahaman,  L.  R.  Dharani,  R.  E.  Mcore 


ONGAMZATION  NAMftS)  AMO  AOO«£SS<ES) 

The  University  of  Missouri-Rolla 
Roila,  MO  65401 


WOIT 


AFOSR-TP. 


omanoahom 


SPONSORNC/ MONITORING  AGENCY  MAME(S)  AMO  AOORESS<ES) 

Dr.  Liselotte  J.  Schioler 

Air  Force  Office  of  Scientific  Research/NE 

Bolling  Air  Force  3ase 

Washington,  DC  20332 


10.  SBONSO0MG/MOMTORMG 


12a.  OtSTMUTKM  /  AVAAAMUTY  STATEMENT 


unlimited 


MDbjMl 


APRi  J  ",ij 


12N.  DISTRIBUTION  COOC 


13.  ABSTRACT  (Mjxumun 200 wotxHl 

In  the  first  reporting  period,  experimental  and  theoretical  studies  were 
performed  to  determine  the  effect  of  key  material  and  processing  parameters  on  the 
sintering  of  model  ceramic  composites. 

Experimentally,  the  effect  of  coarse  inclusions  (ZrC>2)  on  the  sintering  of  a 
fine-grained  matrix  (ZnO)  was  highly  dependent  on  the  forming  method.  Kinetic  data 
and  micros true tural  observations  indicated  that  the  main  factors  controlling 
sintering  were  the  packing  of  the  matrix  immediately  surrounding  the  inclusions  and 
interactions  between  the  inclusions.  A  forming  method  for  alleviating  these  two 
^  ^problems  was  investigated.  Using  this  method,  composites  with  up  to  35  v%  inclusions 
{were  sintered  to  full  density  under  identical  conditions  used  for  the  pure  matrix. 

|  Theoretically,  a  finite  element  program  was  developed  to  describe  the  kinetic 
•data  over  the  entire  sintering  process.  Transient  stresses  were  calculated 
^incrementally  from  the  density- time  data  and  used  to  predict  the  reduction  in  the 
sintering  rate.  The  effects  of  single  and  multiple  inclusions  and  heterogeneity  in 
,  the  matrix  density  were  calculated. 
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1.  STATEMENT  OF  WORK 


Earlier  work  has  shown  that  the  presence  of  dense  second  phase  inclusions  typical  for 
particulate  or  whisker  reinforced  systems  can  cause  substantially  reduced  densification 
rates  in  a  polycrystalline  powder  matrix.  While  a  number  of  explanations  have  been  put 
forward  to  explain  the  reduced  sinterability,  the  mechanism  is  still  ambiguous. 

The  objective  of  the  research  effort  is  to  determine  the  factors  that  control  the 
conventional  pressureless  sintering  of  polycrystalline  ceramic  matrix  composites. 

The  approach  is: 

(i)  invesitgate  how  the  sintering  of  model  composites  is  controlled  by  key  material  and 
processing  parameters  (including  the  volume  fraction,  size  and  distribution  of  the 
reinforcement  phase,  the  uniformity  of  the  matrix,  and  the  sintering  temperature),  and 

(ii)  model  the  observed  behavior. 


2.  STATUS  OF  RESEARCH 


2.1  INTRODUCTION 

This  annual  technical  report  provides  a  description  of  the  work  completed  in  the  first 
year  of  a  three-year  program  whose  aim  is  to  provide  a  basic  understanding  of  the 
processes  that  control  the  sintering  of  polyciystalline  ceram.  natrix  composites  so  that 
materials  with  high  density  and  controlled  microstracture  can  l  .  processed  successfully  by 
conventional  sintering. 

Experimental  studies1-5  have  shown  that  densification  rate  of  a  polycrystalline  powder 
matrix  is  significantly  reduced,  even  at  modest  volume  fractions  of  a  dispersed 
particulate  reinforcement  phase.  The  reduction  in  the  densification  rate  of  glass  matrices 
is  much  smaller  for  an  equivalent  inclusion  content4.  While  a  number  of  explanations5-15 
have  been  put  forward,  the  mechanism  for  the  reduced  sinterability  of  polycrystalline 
matrix  composites  is  still  ambiguous. 

Experimental  and  theoretical  studies  were  carried  out  to  investigate  how  the  sintering 
of  model  composites  is  controlled  by  key  material  and  processing  parameters  (e.g. 
inclusion  volume  fraction,  inclusion  size,  inclusion  distribution,  matrix  uniformity  and 
sintering  temperature).  The  studies  were  initiated  for  a  polycrystalline  composite  system 
consisting  of  a  fine-grained  matrix  (ZnO)  and  an  inert,  rigid,  particulate  inclusion  phase 
(ZrC>2).  ZnO  was  used  as  the  matrix  phase  for  these  initial  studies  because  it  can  be 
sintered  to  full  density  within  the  temperature  range  of  an  existing  dilatometer 
(1600  °C). 

To  provide  further  insight  into  the  processes  influencing  sintering,  green  composite 
bodies  were  formed  by  three  different  processing  techniques:  (a)  conventional  mixing  of 
the  matrix  and  inclusion  powders  followed  by  die-pressing  (hereafter  called  "conventional 
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processing"),  (b)  slip  casting,  and  (c)  coating  the  inclusion  particles  with  the  matrix 
powder  by  in-situ  chemical  precipitation  followed  by  die-pressing  (referred  to  as  "coated 
particles").  For  the  present  report,  it  is  convenient  to  describe  the  experiments  for  the 
three  processing  techniques  separately  (Sections  2.2.1  -  2.23). 

An  important  issue  is  the  effect  of  the  inclusion  phase  on  the  viscosity  of  the 
composite.  The  dependence  of  the  composite  viscosity  on  the  volume  fraction  and  size  of 
the  inclusion  phase  can  provide  insight  into  the  importance  of  inclusion-inclusion  and 
inclusion-matrix  interactions  and  how  thesy  affect  sintering.  It  is  also  important  to 
determine  whether  the  effect  of  an  inclusion  phase  on  the  viscosity  of  a  glass  matrix  is 
qualitatively  different  from  that  for  a  polycrystalline  matrix.  The  effect  of  spherical 
nickel  particles  on  the  viscosity  of  a  soda-lime  glass  matrix  composite  has  been  measured. 
Glass  was  chosen  as  the  matrix  phase  for  the  initial  phase  of  this  study  to  eliminate  the 
complexities  of  grain  growth.  Current  work  is  examining  the  effect  of  ZrCb  inclusions  on 
the  viscosity  of  a  ZnO  matrix  composite.  The  effect  of  the  inclusions  on  the  viscosity  of 
the  glass  matrix  composite  forms  the  fourth  part  of  the  experimental  research  (Section 
2.2.4). 

Theoretically,  modelling  studies  were  initiated  to  account  for  the  experimental  data. 
A  finite  element  program  has  been  developed  to  describe  the  data  over  the  entire  sintering 
process.  The  sintering  problem  is  considered  as  a  special  thermomechanical  problem  in 
which  the  instantaneous  free  strain  rate  as  well  as  the  material  properties  are  density 
dependent.  Experimental  data  for  the  relative  density  versus  time  for  the  composite  is 
used  to  calculate  the  stresses  that  lead  to  a  reduction  in  the  densification  rate. 

However,  the  density  of  the  composite  is  not  known  a  priori.  The  analysis  has  to  be 
performed  in  an  incremental  fashion.  The  density  calculated  from  the  previous  time  step  is 
used  to  calculate  the  instantaneous  stress  and  the  material  parameters  of  the  current  time 
step.  Both  single  inclusion  and  multiple  inclusion  cases  have  been  considered. 
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Inhomogeneity  in  the  matrix  density  has  also  been  considered.  The  modelling  studies  are 
described  in  Section  23. 

The  main  accomplishments  of  the  experimental  and  theoretical  research  are  summarized 
in  Section  2.4. 

The  experimental  and  theoretical  studies  conducted  in  the  first  year  of  the  present 
program  are  described  in  the  following  sections. 
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22  EXPERIMENTAL  STUDIES 


2.2.1  Sintering  of  Model  Ceramic  Composites  Formed  bv  Conventional  Processing 
2.2.1(a)  Experimental  Procedure 

A  fine  grained  ZnO  powder  (Reagent  Grade,  Mallinckrodt  Inc.,  Paris,  KY;  average 
particle  size  ~  0.3  fim)  was  used  as  the  matrix  phase  and  the  inclusion  phase  consisted  of 
relatively  coarse  ZrC>2  particles  (SC  Grade,  Magnesium  Electron  Inc.,  Remington,  NJ)  with 
three  different  average  sizes  (1, 3,  and  14  nm).  The  ZnO  powder  was  the  same  as  that  used 
in  earlier  work  of  Rahaman  and  De  Jonghe.M4,15  Composite  samples  containing  up  to  22.5 
v%  inclusions  (based  on  the  fully  dense  material)  were  formed  by  conventional  processing 
techniques.  The  matrix  powder  and  inclusions  were  mixed  by  milling  in  isopropanol  in  a 
polyethylene  container  using  high  purity  zirconia  balls  as  the  milling  medium.  The  mixture 
was  then  stir-dried  and  ground  lightly  with  a  mortar  and  pestle.  Green  composite  bodies  (6 
mm  in  diameter  by  6  mm)  with  a  matrix  density  of  0.51  ±  0.01  were  formed  by  uniaxial 
pressing  in  a  die  at  pressures  of  »  25  -  30  MPa. 

The  compacted  samples  were  sintered  either  under  isothermal  conditions  at  three  fixed 
temperatures  of  750, 850  and  950  °C,  or  at  a  constant  heating  rate  of  4  °C/min.  Most  of 
the  experiments  were  terminated  after  3  h  for  isothermal  sintering  or  1200  °C  for  constant 
heating  rate  sintering,  but  additional  samples  were  sintered  for  longer  times  (equivalent 
to  higher  temperatures  for  constant  heating  rate  sintering)  to  produce  composites  with 
higher  density.  All  experiments  were  performed  in  air  in  a  dilatometer  (Dilatronic  1600 
°C,  Theta  Industries  Inc.,  Port  Washington,  NY)  that  allowed  continuous  monitoring  of  the 
axial  shrinkage.  The  composite  density  was  measured  from  the  axial  shrinkage  and  the 
initial  density.  The  density  of  the  matrix  phase  of  the  composite  was  calculated  from  the 


6 


measured  composite  density  using  the  method  outlined  earlier.14  The  final  density  of 
selected  sintered  samples  were  also  checked  by  Archimedes  method. 

To  provide  samples  for  grain  size  measurements,  sintering  was  performed  for  the 
appropriate  time  (isothermal  sintering)  or  to  the  appropriate  temperature  (constant 
heating  rate  sintering).  A  different  sample  was  used  for  each  run.  Data  was  obtained  from 
scanning  electron  micrographs  of  polished  and  chemically  etched  surfaces.  Dilute  acetic 
acid  was  used  as  the  etchant.  The  grain  size  was  taken  as  1.5  times  the  mean  intercept 
length  determined  from  about  50  measurements  for  each  sample. 


2.2.1(b)  Results  and  Discussion 


As  outlined  earlier,  sintering  experiments  were  performed  either  at  a  constant  heating 
rate  (4  °C/min)  or  isothermally.  The  results  for  the  constant  heating  rate  experiments 
will  be  considered  first. 

Figure  1  shows  the  data  for  the  relative  density  of  the  composite  matrix,  pm  vs 
temperature  for  ZnO  containing  0-22.5  v%  Zr(>2  inclusions  (average  particle  size  *  14  pm) 
during  sintering  to  1200  °C.  Each  curve  is  the  average  of  two  runs  under  identical 
conditions  and  the  data  are  reproducible  to  within  *  1%.  The  increase  in  the  volume 
fraction  of  the  inclusions,  vj,  caused  a  decrease  in  the  density  at  any  temperature.  The 
composite  density  did  not  reach  a  limiting  value  when  the  unreinforced  ZnO  stopped 
densifying  (at  —  1100  °C);  the  density  continued  to  improve  at  higher  temperature. 

The  matrix  densification  rate,  e  m,  equal  to(l/pm)dpm/dt,  was  calculated  from  the  data 
of  Fig.  1  by  fitting  smooth  curves  to  the  data  and  differentiating;  the  results  are  shown 
in  Fig.  2  as  a  function  of  temperature.  Two  temperature  regimes  are  observed:  one  between 
»«  650  and  950  °C  where  the  densification  rate  at  any  temperature  decreased  with  increasing 
vj  and  the  other  above  *  950  °C  where  the  densification  rate  increased  with  increasing  vj. 
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Fig.  1  .  Matrix  relative  density  versus  temperature  for  conventionally 
processed  ZnO  containing  0-22.5  v%  ZrCb  (average  size  14  microns)  during 
constant  heating  rate  sintering  at  4°C/min.  to  1200°C. 
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Fig.  2  :  Matrix  densification  rate  versus  temperature  for  the  composites 

described  in  Fig.  1. 


Figure  3  shows  the  data  for  both  pm  and  the  matrix  grain  size,  G,  for  ZnO  and  a 
composite  containing  10  v%  Z1O2.  It  is  seen  that  at  any  temperature,  the  inclusions  caused 
a  decrease  in  both  pm  and  G. 

The  data  for  G  vs  pm  is  shown  in  Fig.  4  for  ZnO  and  a  ZnO/10  v%  Zr02  composite.  For 
pmup  to  -  0.92,  the  grain  size  vs  density  relationship  is  approximately  the  same  for  ZnO 
and  the  composite.  This  is  an  interesting  result  and  shows  that  the  grain  size  of  the 
matrix  depends  only  on  the  density  of  the  sample  regardless  of  the  temperature  required  to 
achieve  that  density. 

The  effect  of  inclusion  size  on  the  sintering  of  ZnO  is  shown  in  Fig.  5  for  composites 
containing  10  v%  Zr02  inclusions  with  average  sizes  of  *•  14, 3  and  1  pm.  The  matrix 
density  decreased  with  decreasing  inclusion  size.  This  decrease  cannot  be  explained  by 
current  models. 

For  the  isothermal  sintering  experiments,  Fig.  6  shows  the  data  for  /?m  vs  time  for  ZnO 
and  a  ZnO/10  v%  Z1O2  composite  at  three  different  temperatures:  750, 850  and  950  °C.  The  . 
inclusion  size  was  14  nm.  The  composite  reached  a  limiting  density  which  was  significantly 
lower  than  the  density  of  the  the  unreinforced  sample.  This  behavior  is  qualitatively 
different  from  that  observed  earlier  for  the  constant  heating  rate  experiments  where  the 
matrix  density  increased  with  temperature  and  approached  the  value  for  the  fully  dense 
ZnO. 

The  G  vs  /?m  data  for  isothermal  sintering  at  850  °C  is  shown  in  Fig.  7  for  ZnO  and  a 
ZnO/10  v%  Zr02  composite.  The  data  for  the  composite  is  nearly  identical  to  that  for  the 
unreinforced  ZnO  up  to  -  0.80;  this  value  of  pm  corresponded  to  the  limiting  density 
reached  by  the  composite  matrix.  Further  prolonged  sintering  caused  very  slow  grain  growth 
with  almost  no  increase  in  density. 

To  compare  the  data  from  the  isothermal  and  constant  heating  rate  experiments,  the 
ratio  of  the  densification  rate  of  the  composite  matrix,  e  m,  to  the  densification  rate  of 
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Fig.  3  :  Matrix  relative  density  and  grain  size  versus  temperature  for  ZnO 
and  ZnO/10v%  Zr02  composite  during  constant  heating  rate  sintering  at 
4°C/min.  to  1200°C. 
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Fig.  4  :  Grain  size  versus  matrix  density  for  ZnO  and  ZnO/10v%ZrO2 

composite  during  constant  heating  rate  sintering. 
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Relative  Density 


Fig.  5  :  Matrix  relative  density  versus  temperature  for  ZnO  and  ZnO/10v% 

Zr02  composite  showing  the  effect  of  inclusion  size. 
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Fig.  6  :  Matrix  relative  density  versus  time  for  ZnO  and  ZnO/10v%  ZrCb 

composite  during  isothermal  sintering  at  750,  850,  and  950°C 
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Fig.  7  Grain  size  versus  matrix  density  for  ZnO  and  ZnO/10v%ZrO2 
composite  during  isothermal  sintering  at  850°C. 


15 


the  unreinforced  ZnO,  <  u,  is  shown  in  Fig.  8  as  a  function  of  time  (isothermal  sintering) 
and  as  a  function  of  temperature  (constant  heating  rate  sintering).  The  ratio  €m/€u  was 
almost  constant  with  time  (or  density)  but  increased  with  the  sintering  temperature  for 
isothermal  sintering.  The  constant  heating  rate  sintering  data  increased  with  temperature. 

As  outlined  recently,^  the  data  of  Fig.  8  can  be  expressed  by: 

*m~F<u  (1) 

where  F  is  factor  that  is  strongly  dependent  on  temperature  but  is  weakly  dependent  on  the 
matrix  density.  The  increase  in  F  with  temperature  indicates  that  improved  processing 
conditions  for  ceramic  composites  should  prevail  at  increased  sintering  temperatures.  The 
temperature  dependence  of  F  does  not  follow  from  a  treatment  of  the  viscous  backstress  in 
a  concentric  sphere  model  for  a  densifying  composite.^  Instead  it  indicates  that  an 
additional  process  is  involved  that  is  strongly  dependent  on  temperature. 

Scanning  electron  micrographs  of  (a)  the  unreinforced  ZnO  and  (b)  a  ZnO/lO  v%  Zr02 
composite  which  were  sintered  under  identical  conditions  (4  °C/min  to  1200  °C)  are  shown 
in  Fig.  9.  The  ZnO  is  almost  fully  dense  with  a  nearly  equiaxed  microstructure.  In 
contrast,  the  composite  is  highly  porous  with  the  porosity  situated  in  two  distinct 
regions:  (i)  immediately  surrounding  the  inclusion  particles  and  (ii)  nearly  evenly 
distributed  within  the  bulk  of  the  matrix.  Figure  10  is  a  low  magnification  micrograph  of 
the  composite  in  Fig.  9(b)  and  shows  the  general  distribution  of  the  porosity  immediately 
surrounding  the  inclusions. 

The  sintering  data  and  the  micrographs  described  earlier  indicate  that  the  reduced 
sinterability  of  the  composite  is  due  to  two  main  factors:  one  that  leads  to  highly  porous 
regions  immediately  surrounding  the  inclusion  particles,  and  the  other  that  produces  a 
reduction  in  the  density  (and  grain  size)  in  the  bulk  of  the  matrix  phase.  The  first 
factor  appears  to  be  gross  packing  inhomogeneities  near  the  inclusions  due  to  inefficient 
mixing  between  the  matrix  powder  and  the  inclusion  particles.  The  second  appears  to  be  a 
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Fig.  8  :  Ratio  of  the  matrix  densification  rate  for  the  ZnO/10v%  Zr02 
composite  to  the  densification  rate  for  unreinforced  ZnO  versus  time 
(isothermal  sintering)  or  versus  temperature  (constant  heating  rate 
sintering). 
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(b) 


Fig.  9  :  Scanning  electron  micrographs  of  (a)  ZnO  and  (b)  ZnO/10v%  ZrC>2 

composite,  sintered  under  identical  conditions  (4°C/min.  to  1200°C). 
The  green  bodies  were  formed  by  conventional  processing. 


16 


Fig.  10  :  Low  magnification  scanning  electron  micrograph  of  the  ZnO/10v% 
ZrC>2  composite  of  Fig.  9(b),  showing  regions  of  higher  porosity 

around  the  inclusions. 
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constraint  on  the  matrix  phase  due  to  the  inclusions;  this  constraint  causes  a  reduction 
in  the  sintered  density  (or  an  increase  in  the  porosity)  of  the  composite  matrix  relative 
to  the  unreinforced  matrix.  Since  the  grain  size  of  the  matrix  is  controlled  primarily  by 
the  density,  as  observed  earlier,  therefore  the  lower  matrix  density  leads  to  a  lower 
grain  size.  While  the  exact  nature  of  the  constraint  on  the  matrix  phase  is,  as  yet, 
unknown,  experiments  in  the  present  research  program  indicate  that  it  may  be  due  to 
particle-particle  interactions  between  the  inclusions  (see  Section  2.2.4:  Effect  of  Rigid 
Inclusions  on  the  Sintering  and  Viscosity  of  Glass). 

In  summary,  the  work  outlined  so  far  indicates  that  improved  sintering  of  ceramic 
composites  should  be  achieved  with: 

(i)  the  use  of  higher  sintering  temperatures,  and 

(ii)  processing  techniques  which  maximize  the  packing  of  the  matrix  powder  around  the 
inclusion  phase,  and  which  reduce  the  constraint  on  the  matrix  phase  due  to  the  inclusion 
particles. 
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2.2.2  Sintering  of  Model  Ceramic  Composites  Formed  bv  Slip  Casting 


2.2.2(a)  Experimental  Procedure 

The  matrix  powder,  ZnO,  and  the  inclusion  particles,  Z1O2,  were  the  same  as  those  used 
earlier  for  the  composites  formed  by  conventional  processing  (i.e.  Section  2.2.1).  The  ZnO 
powder  and  the  Z1O2  inclusions  (average  particle  size  -  14  or  3  /im)  were  dispersed  in 
water  at  a  pH  of  8.5  to  produce  a  slurry  with  a  solids  content  of  -  80  w%.  Damn  811  was 
used  as  a  dispersant.  The  slurry  was  cast  into  gypsum  molds  to  produce  samples  (6  mm  in 
diameter  by  6  mm)  for  sintering.  The  sintering  experiments  and  microstructural 
observations  were  performed  in  the  same  way  as  described  earlier  for  the  conventionally 
processed  composites. 

2.2.2(b)  Results  and  Discussion 

Figure  11  shows  the  data  for  the  matrix  relative  density  vs  temperature  for  the 
unreinforced  ZnO  and  a  ZnO/10  v%  Zr02  composite  during  constant  heating  rate  sintering  at 
4  °C/min  to  1300  °C.  Compared  to  the  conventionally  processed  bodies,  the  unreinforced  ZnO 
started  densifying  at  a  somewhat  higher  temperature;  this  difference  in  the  sintering 
behavior  may  be  due  to  a  small  amount  of  hydrolysis  of  the  ZnO  during  slip  casting.  Since 
the  green  bodies  were  formed  by  a  nearly  identical  process,  therefore  the  effect  of  this 
hydrolysis  would  be  expected  to  affect  the  unreinforced  matrix  and  the  composite  in 
approximately  the  same  way.  The  unreinforced  ZnO  body  eventually  reached  nearly  full 
density  at  *  1250  °C.  It  is  seen  from  Fig.  11  that  the  matrix  density  of  the  composite  at 
any  temperature  is  much  lower  that  for  the  unreinforced  ZnO.  The  highly  reduced  density  of 
the  composite  matrix  is  comparable  to  the  reduction  observed  for  the  conventionally 
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Fig.  11  :  Matrix  relative  density  for  slip  cast  ZnO  and  ZnO/10v%  Zr02 
composite  during  constant  heating  rate  sintering  at  4°C/min.  to  1300°C. 
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processed  samples. 

The  ratio  of  the  matrix  densification  rate  of  the  composite  to  the  densification  rate 
of  the  unreinforced  ZnO  (i.e.  e  m/«  u)  as  a  function  of  temperature  is  shown  in  Fig.  12.  The 
curve  has  almost  the  same  shape  but  is  shifted  to  somewhat  higher  temperatures  compared  to 
the  results  for  the  conventionally  processed  bodies  (Fig.  8). 

Scanning  electron  micrographs  of  (a)  the  unreinforced  ZnO  and  (b)  a  ZnO/lO  v%  Zr02 
composite  which  were  sintered  under  identical  conditions  (4  °C/min  to  1300  °C)  are  shown 
in  Fig.  13.  It  is  seen  that  (i)  the  composite  matrix  is  more  porous  than  the  unreinforced 
sample  and  (ii)  the  packing  of  the  matrix  phase  close  to  the  inclusions  is  much  better 
compared  to  the  conventionally  processed  composites  (Fig.  9). 

The  kinetic  data  indicate  that  the  sinterability  of  the  slip  cast  composites  is  not 
significantly  improved  compared  to  the  conventionally  processed  composites.  However,  the 
packing  uniformity  of  the  matrix  phase  immediately  surrounding  the  inclusion  particles  is 
much  improved.  It  is  worth  comparing  the  green  composite  microstructure  expected  from  the 
conventional  processing  and  slip  casting  routes.  Both  conventional  processing  and  slip 
casting  are  expected  to  lead  to  a  random  distribution  of  the  inclusion  particles  within 
the  matrix.  Slip  casting  is,  however,  expected  to  produce  better  packing  uniformity  of  the 
matrix  phase  close  to  the  inclusions  and  in  the  rest  of  the  matrix.  The  feature  that  is 
common  to  both  processing  techniques,  therefore,  is  the  random  distribution  of  the 
inclusion  particles  within  the  matrix.  Hence,  for  the  slip  cast  bodies,  particle-particle 
interactions  that  result  from  the  random  distribution  of  the  inclusion  phase  might  be  the 
main  factor  causing  the  reduced  sinterability  of  the  slip  cast  samples.  For  conventionally 
processed  composites,  an  additional  factor  is  the  packing  uniformity  around  the  inclusion 
phase. 
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Fig.  12 :  Ratio  of  the  matrix  densification  rate  for  the  slip  cast  ZnO/10v% 
ZK)2  composite  to  the  densification  rate  of  the  unreinforced  ZnO  versus 
temperature. 
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(b) 

Fig.  13  :  Scanning  electron  micrographs  of  (a)  ZnO  and  (b)  ZnO/10v%  ZrC>2 

composite,  sintered  under  identical  conditions .  The  green  bodies 
were  formed  by  slip  casting. 
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As  outlined  earlier,  the  data  for  the  conventionally  processed  composites  and  slip 
cast  composites  indicate  that  the  matrix  powder  padring  immediately  surrounding  the 
inclusions  and  partide-partide  interactions  between  the  inclusion  partides  (which  are 
randomly  distributed  in  the  matrix)  might  be  the  important  factors  responsible  for  the 
reduced  sinterabilty  of  ceramic  particulate  composites.  The  objective  of  this  phase  of  the 
experimental  work  was  to  investigate  whether  these  problems  can  be  alleviated  by  improved 
processing.  The  forming  method  described  in  this  section  involved  coating  the  individual 
inclusion  particles  with  a  cladding  of  the  matrix  powder  by  chemical  precipitation, 
followed  by  die-pressing  of  the  coated  particles. 

2.23(a)  Experimental  Procedure 

The  Z1O2  inclusion  particles  (average  size  14, 3  and  1  /im)  were  the  same  as  those  used 
earlier  in  the  formation  of  composites  by  conventional  processing  and  by  slip  casting.  The 
inclusion  particles  were  coated  with  a  matrix  phase  of  ZnO  by  chemical  preripitation.  In 
the  coating  process,  the  Z1O2  particles  were  dispersed  in  a  solution  of  zinc  chloride  and 
urea  (Aldich  Chemical  Company)  at  a  pH  of  6.7.  All  solutions  were  prepared  in  distilled 
water.  The  system  was  aged  at  70  °C  with  continuous  stirring  for  *  20  hours,  cooled  to 
room  temperature  and  the  solids  were  recovered  by  filtration.  After  washing  with  distilled 
water,  the  solids  were  dried  at  150  °C  for  24  hours.  Coaled  particles  containing  up  to  =* 

40  v%  inclusion  particles  (based  on  the  fully  dense  composite)  were  synthesized.  For 
comparison,  pure  ZnO  powder  (i.e.  containing  no  inclusion  particles)  was  synthesized  by 
the  same  chemical  prcipitation  process. 

The  powders  were  compacted  by  uniaxial  pressing  in  a  die  to  produce  green  bodies  with 


approximately  the  same  matrix  density  (0.52  -  0.02).  The  compacts  were  sintered  in  air  at 
a  constant  heating  rate  of  4  °C/min  to  1500  °C  in  the  dilatometer  used  earlier. 

2 23  Results  and  Discussion 

Figure  14  shows  scanning  electron  micrographs  of  (a)  ZrCh  inclusion  particles  with  an 
average  size  of  *  3  pm,  and  (b)  the  same  particles  after  they  were  coated  with  ZnO  powder 
by  chemical  precipitation.  The  coated  particles  appear  somewhat  agglomerated.  X-ray 
diffraction  and  energy  dispersive  X-ray  analysis  in  the  scanning  electron  microscope  for 
coated  particles  containing  20  v%  of  the  iclusion  phase  did  not  reveal  the  presence  of  the 
ZrC>2  inclusions.  This  is  a  strong  indication  that  the  inclusion  particles  were  indeed 
coated  by  the  ZnO  matrix  phase.  As  a  further  check,  high  voltage  transmission  electron 
microscopy  will  be  used  in  future  work. 

The  matrix  relative  density  versus  temperature  for  the  composites  formed  from  coated 
particles  is  shown  in  Fig.  15.  The  Z1O2  inclusion  size  was  ~  3  pm.  For  comparison,  the 
data  for  the  ZnO  precipitated  under  the  same  conditions  as  the  coated  particles  is  also 
shown.  The  pure  ZnO  started  to  density  at  a  higher  temperature  compared  to  the  commercial 
ZnO  powder  (particle  size  *  0.3  pm)  used  for  the  conventionally  processed  samples. 
Although  the  particle  size  of  the  chemically  precipitated  powder  has  not,  as  yet,  been 
measured,  the  difference  in  sintering  characteristics  between  the  two  ZnO  powders  might  be 
related  to  a  difference  in  particle  size.  The  data  of  Fig.  15  show  that  composites  with  up 
to  35  v%  of  inclusions  densified  to  nearly  full  density  under  identical  conditions  used 
for  the  unreinforced  matrix  which  reached  a  density  of  only  0.93.  Further,  the 
densification  of  the  matrix  is  rather  insensitive  to  the  inclusion  volume  fraction  for  the 
range  studied. 

Figure  16  shows  data  for  the  matrix  density  as  a  function  of  time  for  composites  with 
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Fig.  15  :  Matrix  relative  density  versus  temperature  for  ZnO  and  ZnO 

containing  up  to  35  V%  Zr02  inclusions  (  particle  size  3  micron) 
during  constant  heating  rate  sintering.  The  composites  were 
formed  from  coated  particles. 
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a  fixed  inclusion  content  (20  v%)  but  with  different  inclusion  sizes  of  14, 3  and  1  ^m.  It 
is  seen  that  the  inclusion  size  has  almost  no  effect  on  the  densification  of  the  matrix. 

The  insensitivity  of  the  matrix  densification  to  inclusion  content  and  inclusion  size 
(Figs.  15  and  16)  is  remarkably  different  from  the  data  for  the  conventionally  processed 
bodies  and  for  the  slip  cast  bodies.  The  results  show  very  clearly  that  the  effects  of 
rigid  inclusions  on  the  sintering  of  ceramic  composites  are  processing-related. 

Figure  17  shows  scanning  electron  micrographs  of  (a)  a  polished  surface  and  (b)  a 
fractured  surface  of  a  composite  containing  20  v%  inclusions  (particle  size  3  n m)  that  was 
formed  from  coated  particles  and  sintered  at  4  °C/min  to  1500  °C.  It  is  seen  that  the 
matrix  (including  the  region  immediately  surrounding  the  inclusions)  is  nearly  fully 
dense.  The  fractured  surface  shows  that  the  matrix  consists  of  nearly  equiaxed  grains. 

In  summary,  the  results  for  the  composites  formed  from  coated  powders  show  that  the 
problems  associated  with  the  sintering  of  ceramic  composites  can  be  alleviated  by  the 
present  technique  of  coating  the  inclusion  particles  with  the  matrix.  Since  the  main 
advantages  of  this  method  are  (i)  the  improvement  of  the  matrix  packing  immediately 
surrounding  the  inclusions  and  (ii)  reduction  in  the  particle-particle  interactions 
between  the  inclusions,  therefore  the  data  provide  further  evidence  that  non-uniform 
packing  of  the  matrix  immediately  surrounding  the  inclusions  and  inclusion  particle 
interactions  are  the  main  factors  responsible  for  the  reduced  densification  rates  in 
conventionally  processed  composites.  Slip  casting  leads  to  an  improvement  in  the  matrix 
packing  immediately  surrounding  the  inclusions  (compared  to  conventional  processing)  but 
the  interactions  between  the  inclusion  particles  are  not  significantly  reduced. 
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Fig.  16  :  Matrix  relative  density  versus  temperature  for  ZnO/20  V%  Zr02 
composites  formed  from  coated  particles,  showing  almost  no 
effect  of  the  Zr02  inclusion  size. 
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Fig.  17  :  Scanning  electron  micrographs  of  (a)  a  polished  surface  and  (b)  a 
fractured  surface  of  a  ZnO/20  V%  Zr02  composite  formed  from 
coated  particles  and  sintered  at  4°C  /  min.  to  1485°C. 
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2.2.4  Effect  of  Rigid  Inclusions  on  the  Composite  Viscosity 


The  experiments  described  earlier  indicated  that  interactions  between  the  inclusion 
particles  provide  an  important  barrier  to  the  sintering  of  composites  formed  by 
conventionally  processing  and  by  slip  casting.  Significant  interactions  between  the 
inclusion  particles  should  lead  to  an  increase  in  the  viscosity  of  the  composite.  An 
important  issue  therefore  is  the  effect  of  the  inclusion  phase  on  the  viscosity  of  the 
composite.  The  dependence  of  the  composite  viscosity  on  the  volume  fraction  and  size  of 
the  inclusion  phase  can  provide  insight  into  the  importance  of  inclusion-inclusion  and 
inclusion-matrix  interactions  and  how  these  affect  sintering.  Since  the  sintering  of  glass 
matrix  composites  is  qualitatively  different  from  that  for  polycrystalline  composites,  it 
is  also  important  to  determine  whether  the  effect  of  an  inclusion  phase  on  the  viscosity 
of  a  glass  matrix  is  different  from  that  for  a  polycrystalline  matrix.  In  the  first  phase 
of  these  experiments,  the  effect  of  spherical  nickel  particles  on  the  viscosity  of  a 
soda-lime  glass  matrix  composite  was  measured.  Glass  was  chosen  as  the  matrix  phase  for 
the  initial  phase  of  this  study  to  eliminate  the  complexities  of  grain  growth.  Current 
work  is  examining  the  effect  of  ZrC>2  inclusions  on  the  viscosity  of  a  ZnO  matrix  composite 
and  the  results  will  be  reported  later. 

2.2.4(a)  Experimental  Procedure 

Spherical  soda-lime  glass  powder  (diameter  10-20  jam;  MOSCI  Corporation,  Rolla,  MO)  was 
used  as  the  matrix  phase  and  the  inclusion  phase  consisted  of  nearly  spherical  Ni 
particles  with  three  different  size  ranges:  6-10, 100-150,  and  150-300  fim  (Johnson 
Matthey,  Ward  Hill,  MA).  The  glass  powder  and  inclusion  particles  were  mixed  in  methanol 
and  stir-dried.  Approximately  5  v%  of  carbowax  was  used  as  a  binder.  For  comparison,  the 
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pure  glass  powder  was  subjected  to  the  same  mixing  treatment.  The  powders  were  then 
pressed  uniaxially  in  a  die  to  form  samples  (6  mm  in  diameter  by  6  mm)  with  a  matrix 
relative  density  of  *  0.60.  The  compacted  samples  were  heated  slowly  to  **  400  °C  to  remove 
the  binder  and  then  presintered  in  argon  gas  at  4  °C/min  to  710  °C  to  produce  compacts 
with  sufficient  green  strength  for  manipulation  in  the  dilatometer.  All  samples  were  given 
the  same  binder  burnout  and  presintering  treatment.  The  presintering  treatment  resulted  in 
an  increase  of  2-3%  in  the  relative  density  of  the  sample. 

The  unreinforced  glass  and  composite  samples  were  sintered  isothermally  with  or 
without  load  in  argon  gas  at  730  °C  in  a  loading  dilatometer^,  17  that  allowed  continuous 
monitoring  of  the  axial  shrinkage.  For  sintering  with  load,  the  initial  stress  on  the 
compact  was  18  kPa.  The  density,  creep  strain  and  creep  viscosity  during  sintering  were 
determined  from  the  data  using  a  procedure  outlined  earlier.  16*17 

In  a  separate  set  of  experiments,  the  creep  behavior  of  nearly  fully  dense  glass  and 
composite  samples  were  also  measured.  The  compacts  were  first  sintered  freely  (i.e. 
without  load)  at  730  °C  to  nearly  full  density.  In  the  creep  experiments,  a  sample  was 
placed  in  the  loading  dilatometer  and  after  the  desired  temperature  (730  °C)  had  been 
reached,  a  stress  of  18  kPa  was  applied.  The  axial  strain  of  the  sample  was  measured  as  a 
function  of  time. 

Microstructural  observations  of  the  samples  various  stages  of  sintering  were  performed 
by  scanning  electron  micoscopy. 

2.2.4(b)  Results  and  Discussion 

Figure  18  shows  the  data  for  the  matrix  relative  density  versus  time  for  the 
unreinforced  glass  and  glass  matrix  composites  containing  10  v%  Ni  inclusions  during 
sintering  at  730  °C.  (The  load  on  the  samples  was  «•  3  kPa,  corresponding  to  the  load 
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Fig.  18  :  Matrix  relative  density  versus  time  for  glass  and  glass  matrix 

composites  containing  10  V%  Ni  particles  sintered  at  730°C  under 
a  uniaxial  load  of  3  kPa. 
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.  required  to  maintain  contact  of  the  dilatometer  pushrod  on  the  sample).  These  density  data 
were  calculated  from  the  measured  axial  strain  data  using  the  procedure  outlined 
earlier.16,17  The  size  of  the  the  Ni  inclusions  was  6-10  pm  in  one  set  of  experiments  and 
150-300  pm  in  the  other  set.  The  data  for  a  third  set  of  composites  with  Ni  inclusion  size 
of  100-150  pm  have  been  omitted  to  maintain  clarity.  It  is  seen  that  both  the  glass  sample 
and  the  composite  with  the  Ni  inclusion  size  of  150-300  pm  reached  nearly  full  density 
within  approximately  the  same  sintering  time.  The  composite  with  the  Ni  inclusion  size  of 
6-10  pm,  however,  reached  a  density  of  -  0.93. 

As  outlined  earlier,  the  creep  of  nearly  fully  dense  samples  was  also  studied.  Figure 
19  shows  the  axial  strain  versus  time  for  the  dense  unreinforced  glass  samples  and 
composites  containing  10  v%  inclusions  with  three  different  sizes.  The  applied  uniaxial 
stress  was  18  kPa.  The  strain  at  any  time  was  smaller  for  the  composite  compared  to  the 
unreinforced  sample  and  decreased  with  decreasing  inclusion  size.  Corresponding  data  for 
.  composites  containing  20  v%  inclusions  with  two  different  sizes  (150-300  and  6-10  pm)  are 
shown  in  Fig.  20.  For  these  composites,  the  sample  with  the  larger  inclusion  size  was 
nearly  fully  dense  but  that  with  the  lower  inclusion  size  had  a  matrix  density  of  only  *• 

0.85. 

Figure  21  shows  the  creep  viscosity  versus  the  matrix  density,  pm,  for  the 
unreinforced  glass  and  the  glass  matrix  composites  containing  10  v%  Ni  inclusions.  (The 
units  for  the  viscosity  are  MPa.min,  and  not  MPa/min  as  printed  on  the  figures.)  Over  most 
of  the  density  range  (pm  up  to  «•  0.90),  the  viscosity  of  the  composite  is  insensitive  to 
the  presence  of  the  inclusions.  Above  pm  -  0.90,  the  viscosity  shows  a  small  dependence  on 
the  inclusion  size  and  this  is  confirmed  by  the  data  for  the  fully  dense  samples.  The 
viscosity  of  the  fully  dense  samples  increases  by  a  factor  of  *» 5  from  the  unreinforced 
glafcs  to  the  composite  with  Ni  inclusion  particles  of  6-10  pm. 

The  creep  viscosity  as  a  function  of  pm  for  the  composites  with  20  v%  Ni  particles  is 
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Fig.  19  :  Axial  strain  versus  time  for  nearly  fully  dense  samples  of  glass 
and  glass  /  10  V%  Ni  composite  subjected  to  a  uniaxial  stress  of 
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Fig.20  :  Axial  strain  versus  time  for  dense  glass  /  20  V%  Ni  composites 
subjected  to  a  uniaxial  stress  of  18  kPa. 
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Fig.21  :  Creep  viscosity  versus  matrix  density  for  glass  and  glass  /  10  V% 
Ni  composites  at  730°C. 


shown  in  Fig.  22.  For  comparison,  the  viscosity  of  the  unreinforced  glass  is  also  shown.  A 
much  stronger  inclusion  size  effect  is  observed.  At  a  density  of  *  0.85,  the  viscosity  of 
the  composite  with  Ni  inclusion  size  of  6-10  iim  is  more  than  2  orders  of  magnitude  greater 
than  the  viscosity  of  the  unreinforced  glass. 

Scanning  electron  micrographs  of  a  glass  composite  containing  10  v%  Ni  inclusions 
(150-300  /xm)  after  (a)  the  pre-sintering  step  and  (b)  sintering  for  42  min  are  shown  in 
Fig.  23.  The  inclusions  are  well  dispersed  within  the  matrix  and  after  sintering,  the 
glass  is  fully  densified  around  the  inclusions.  Corresponding  micrographs  for  the 
composite  with  inclusions  6-10  /xm  in  size  are  shown  in  Fig.  24.  The  inclusions  are 
somewhat  agglomerated  and  after  sintering,  the  glass  was  unable  to  fill  the  pores  within 
the  agglomerates. 

The  data  show  that  for  the  inclusion  size  of  6-10  /xm  (which  is  comparable  to  the 
particle  size  of  the  glass  powder)  the  viscosity  of  the  composite  increases  dramatically 
in  the  range  of  10-20  v%  inclusions.  These  results  are  somwhat  similar  in  nature  to  the 
theoretical  predictions  of  Scherer  and  Jagota^.  They  indicate  that  inclusion  particle 
interactions  are  responsible  for  the  reduced  densification  rates  above  10  v%.  It  would  be 
interesting  to  determine  how  the  viscosity  of  polycrystalline  composites  depends  on 
inclusion  size  and  volume  fraction. 
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Fig.2?  :  Creep  viscosity  versus  matrix  density  for  glass  an  glass  /  20  V% 
Ni  composites,  compared  to  the  viscosity  for  the  unreinforced 
glass. 
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Fig.  24  :  Scanning  electron  micrographs  of  a  glass  / 10  V%  Ni  composite 
(Ni  particle  size  6-10  micron)  after  (a)  pre-firing  for  20  min  at 
710°C  and  (b)  after  sintering  42  min  at  730°C. 


2.3  THEORETICAL  MODELING 


2.3. 1  Background. 

It  has  been  observed  that  the  presence  of  coarse,  inert,  rigid  inclusions  in  a  fine 
'  grained  matrix  drastically  reduces  the  densiflcation  rates  in  the  sintering  process  [1,3] . 

A  significant  theoretical  modeling  work  has  been  done  in  trying  to  explain  this 
phenomena.  The  reasons  have  been  generally  attributed  to  the  development  of  transient 
stresses  and  the  sintering  matrix  is  commonly  assumed  to  be  viscoelastic.  A  commonly 
used  model  system  consists  of  rigid  inclusions  and  viscoelastic  matrix,  and  viscoelastic 
backstresses  are  predicted  to  cancel  part  of  the  sintering  stress  (sintering  potential)  so 
that  the  matrix  densification  rate  is  reduced  [6,7]  .  However,  Scherer's  theory  [5]  and 
De  Jonghe  and  Rahaman's  analysis  [13]  showed  that  the  backstresses  were  too  small  to 
account  for  the  drastic  reduction  in  densification  rate. 

A  number  of  suggestions  have  been  put  forward  to  account  for  the  drastic  reduction 
in  the  densification  rate.  Lange  [9]  proposed  network  formation  of  inclusions.  However 
network  formation  is  not  significant  at  inclusion  volume  fraction  below  15%  while 
experiments  showed  that  the  densification  rates  of  polycrystalline  composite  were 
drastically  reduced  even  at  low  inclusion  volume  fraction  (<  15%).  Lange  [8,21]  also 
attributed  reduction  in  densification  rate  to  damage,  e.g.  voids,  in  the  matrix  produced  by 
nonuniform  densification.  At  present,  no  detail  analysis  has  been  done  for  the  damage 
formation  ,  development  and  its  effect  on  the  reduction  in  densification  rates. 

A  number  of  other  factors  that  may  lead  to  the  drastic  reduction  in  densification 
rates  in  polycrystalline  matrix  composites  suggested  by  Scherer  and  Bordia  [11].  These 
include  the  competition  between  densification  and  coarsening  ,  and  the  development  of 
anisotropies  resulting  from  the  unequal  stresses  in  radial  and  tangential  directions 
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experienced  by  the  shrinking  matrix.  Very  recently,  some  efforts  have  been  made  in 
analysing  each  of  these  problems  by  using  both  analytical  and  numerical  methods.  Jagota 
developed  a  finite  element  analysis  for  neck  growth  between  two  viscous  spherical 
particles  [28]  and  also  used  truss  structures  in  simulating  the  contacts  among  ceramic 
powder  particles  [25]  .  A  semi-analytical  solution  of  a  viscous  media  between  two  rigid 
spheres  is  given  by  Andareson  [24]  .  But  no  single  factor  has  been  found  to  be  significant 
alone. 


2.3.2  Necessity  for  developing  a  finite  element  program. 

The  objective  of  the  present  research  is  to  try  to  approach  some  problems  for  which 

the  closed  form  solutions  are  mathematically  untractable.  A  finite  element  program  is 

developed  to  describe  the  constrained  sintering  matrix  through  the  entire  sintering 

process.  Although  the  heating  schedule  for  a  sintering  process  can  be  quiet  different, 

isothermal  sintering  has  drawn  much  attention  since  temperature  effect  can  be  excluded. 

The  special  situation  of  the  isothermal  sintering  problem  is  that  the  instantaneous  free 

strain  rate,  as  well  as  the  material  properties,  can  be  treated  as  density  dependent.  The 

relation  between  relative  density  and  time  derived  from  experimental  data  of 

unconstrained  sintering  is  used  to  calculate  the  sintering  strain  loads  which  in  turn  is  used 

in  predicting  the  reduction  in  densification  rate  of  constrained  sintering  matrix  as  shown 

in  Fig.25.  The  sintering  strain  load  which  can  be  treated  as  a  kind  of  sintering  potential 

dp 

in  equivalence  can  be  expressed  in  terms  of  relative  density  p  and  its  time  derivative 
of  the  unconstrained  sintering.  That  is 


1  dp 
&f~  3  p  dt 


(2.3.1) 


However  the  density  of  the  contrained  sintering  matrix  during  the  sintering  process  is  not 
known  apriori.  This  kind  of  nonlinear  problem  has  to  be  approached  in  an  incremental 


45 


0.8- 


•  • 

Ajjsuaa  9Aj)B|ay 


46 


Experimental  data  of  relative 
density  vs  time. 


fashion.  The  density  calculated  from  the  previous  increment  is  used  in  determining  the 
instantaneous  sintering  strain  load,  as  well  as  the  material  properties,  in  the  current 
increment,  under  the  assumption  that  both  the  contrained  and  unconstrained  matrices 
have  the  same  sintering  potential  at  the  same  density. 

Experience  with  NASTRAN.  During  the  early  part  this  research  project,  the 
authors  used  MSC/NASTRAN  program's  viscoelastic  analysis  function  in  attempting  to 
deal  with  the  sintering  problem.  NASTRAN  is  a  large  finite  element  analysis  package 
program  and  has  a  viscoelastic  analysis  function.  But  our  experience  showed  that  it  is 
not  suitable  for  simulating  the  sintering  process.  The  main  problem  is  that  the  thermal 
load  term  which  is  treated  as  sintering  strain  load  by  taking  the  thermal  coefficient  as  a 
linear  free  strain  rate  and  the  temperature  variable  as  time.  The  incompatibility  with  the 
sintering  model  is  that  the  free  strain  load  can  not  be  linked  to  density  but  merely  time. 
Referring  to  Fig.  2S,  for  example,  at  time  /(  the  relative  density  of  the  uncontrained 
matrix  is  p„  at  point  A  while  the  relative  density  of  the  constrained  matrix  is  pe  at  point 
C.  The  slope  p  used  in  equation  (2.3.1)  for  calculating  the  free  strain  load  should  be  the 
one  at  point  B  corresponding  to  time  u  not  at  A  corresponding  to  time  This  shows  how 
the  density  dependent  free  strain  rate  differs  from  the  time  dependent  one.  Since  the 
calculated  density  from  the  previous  step  can  not  be  automatically  fed  back  for  use  in  the 
next  step  in  NASTRAN,  the  density  dependent  sintering  strain  load  can  not  be  satisfied, 
nor  the  density  dependent  material  property.  These  limitations  of  NASTRAN 
necessitated  the  development  of  a  new  finite  element  program  for  the  propose  of 
simulating  sintering  process. 

2.3.3  Viscoelastic  Formulation. 


4? 


General  Model.  In  viscoelasticity  theory,  there  are  two  general  models  [19,32]  . 


The  general  Maxwell  model  as  shown  in  Fig.26,  consists  of  multiple  simple  Maxwell 
elements  in  parallel.  We  consider  one  of  them  for  studying  the  sintering  behavior 
Its  total  strain  rate  is 

{*}  =  {e‘}  +  {fiv}  +  {e*}  (2.3.2) 


where 

(e)  total  strain  rate 

(s')  elastic  strain  rate 

{s’}  viscous  strain  rate 

(e’}  sintering  strain  rate 


For  isotropic  viscoelastic  material,  the  viscous  component  of  strain  rate  {s’}  can  be  divided 
into  deviatoric  part  {£•'}  and  volumetric  part  {e*}  .  As  in  isotropic  elastic  case 


{O  =■  { i “)  +  (€*}  =  -±-  {»„}  +  !Jm}  (2.3.3) 

-In 

where  and  tjk  are  shear  and  bulk  viscosities,  respectively.  So  that  in  3  dimension, 
equation  (2.3.2)  can  be  expressed  as 


where  a(  =  —  p/3p)  is  linear  free  strain  rate.  Let  shear  and  bulk  elastic  moduli  be  denoted 
as  ft  and  K  respectively,  then  shear  and  bulk  relaxation  times  rM  and  xK  are  defined  as: 
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Fig. 26  General  Maxwel  M  odel 


ft  ’ 


_  *Ik 

T*~T 


(2.3.5) 


and  the  deviatoric  and  volumetric  stresses  are 
ai  =  (2*„  -<ry-  <rs)/3, 
ai  =  (2 (Ty  -<r.  -  <r*)/3, 

^  =  (2<r*  -  Ok  -  *,)/3, 

<fm  -  («r,  +  <ry  +  «t£)/3 

After  some  matrix  manipulation  with  elastic  stiffness  matrix,  the  following  relation  is 
obtained 

W-[0](i}-3*:(i‘)-^(«')-^W  .  (2.3.6) 

where  [Z>]  is  the  elastic  stiffness  matrices. 

For  the  general  Maxwell  model  consisting  of  multiple  simple  Maxwell  elements  in 
parallel,  the  expression  becomes  a  summation  over  each  of  these  simple  Maxwell 
elements 


W  -  (£[»,])«  -  ^(3*$}  +  -fL  (4  +  _L  (^)) 


(2.3.7) 


Similar  expressions  for  the  general  Kelvin  model,  shown  in  Fig.27,  consisting  of  one 
simple  Maxwell  element  and  multiple  simple  Kelvin  elements  in  a  series  can  also  be 
obtained.  However,  for  most  materials,  it  is  very  difficult  to  obtain  parameters  for  the 
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multiple  element  models.  An  efficient  experimental  technique  of  obtaining  the  material 
parameters  even  in  the  three  or  four  element  standard  solid  model  for  ceramics  is  yet  to 
be  developed  [30]  .  It  is  therefore  practical  to  incorporate  only  the  standard  solid  model 
into  the  finite  element  formulation  as  done  in  several  package  programs,  such  as 
NASTRAN  and  ADINA. 


Standard  Solid  Model.  The  standard  solid  model  consists  of  one  simple  Maxwell 
element  and  one  simple  Kelvin  element  as  shown  in  Fig.  28.  For  a  Kelvin  element, 


(»')  =  W*' )  +  <»„>  - 


(2.3.8) 


hence 


(2.3.9) 


In  the  standard  solid  model, 


{£*}  =  -  J-  {«“}  +  (  V-  + -T-  V}  -  ~  {£rf} 

2/i  V  2*m  J  * 


3 ij  {<)  +  (  3J*  +  3^  ){</}  -  i  {£"} 


(2.3.10) 


(2.3.11) 


tf)  -  Wrt  -  /<(  i  +  ^7  )('')  +  2,(  -£■  (a"}) 


(2.3.12) 


!<)  “  3*{iJ  -  *(  TJ7  +  )  W  +  3«(  -fi-  {:„})  +  3A'{i'}  (2.3.13) 


simpler  models  can  be  obtained  by  setting  any  of  /i,  K,  rj„  tjK,  n*  A':,  tjM  and  qKl  to  be 
either  infinity  or  zero. 
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2.3.4  Finite  Element  Program.. 


This  computer  code  is  used  for  viscoelastic  analysis  of  the  sintering  matrix  with 
inclusions  by  the  finite  element  method.  A  brief  description  of  a  simple  Maxwell  model 
in  two  dimension  (plane  strain)  is  given  below. 

Solution  Procedure.  At  the  beginning  of  nth  time  step  ,  the  volume  around  each 
sampling  point  (usually  Gaussian's  point)  is  calculated  by  the  Gaussian's  numerical 
integration 


^  =  J(c,C)w 


(2.3.14) 


where  the  J(c,  C)  is  the  Jacobian's  matrix  mapping  the  element  from  the  global 
coordinates  to  the  local  coordinates,  the  w.  and  kv  are  the  weight  coefficients  at  the 
sampling  point  in  local  coordinates  (q,  Q  • 

Then  the  relative  density  averaged  over  that  sampling  volume  can  be  calculated  as 


ft 

P 


(2.3.15) 


The  linear  free  strain  rate  is  calculated  as 


1  dpu{p) 

Sf~  3  p  dt 

where  pu  is  the  relative  density  of  the  unconstrained  matrix. 
The  constitutive  parameters  are  all  density  dependent, 


(2.3.16) 


E"'=  E(pn)  (2.3.17) 

vn  =  v(pn)  (2.3.18) 
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nn=n(pn) 


(2.3.19) 


then  the  shear  modulus  and  relaxation  time  can  be  derived  as 


n 

t*  = 


E" 

2(1  +  v") 


n 

T  = 


It 

P 


The  current  coordinates  of  nodes  are 


{Xn}  =  {y1-1}  +  {Ay-1} 

The  current  strain-displacement  relation  matrix  is 


[Bn]  =  [5({x"})] 


The  current  elastic  stress-strain  relation  matrix  is 


=  [/>(/,  /)] 


For  convenience,  a  coefficient  of  the  simple  Maxwell  model  is  denoted  as 


kn  = 


At 

,+T7 


The  equation  of  equilibrium  [29]  is 


+/  =  0 
n 


(2.3.20) 

(2.3.21) 

(2.3.22) 

(2.3.23) 

(2.3.24) 

(2.3.25) 

(2.3.26) 
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The  incremental  form  of  equilibrium  equation  is 


f  lBnf{±an}in  +  \f  =  0 

Jo 


(2.3.27) 


The  viscoelastic  stress-strain  relation  of  simple  Maxwell  model  in  plane  strain  is 


{A«"}  =  rcCM/THA*"}  -  (1  +  v"){A£;»  -  ( )*  V} 


(2.3.28) 


where 


{</■}  =  {ct""1}  +  {Act"-'} 


(2.3.29) 


Substituting  the  {Act"}  into  the  incremental  form  of  equilibrium  equation  (2.3.27)  yields 


W  [5fl]r[Z>',][£',]</n]  {AW"}  =  r(l +  /’)('  lBnf[Dn-]{±s}}dQ. 

Jo  /  Jo 


r  Jo 


(2.3.30) 


where  <l>"  is  called  the  residual  force 


* 

<t>"  =  [B^/i^dn+fto 

Jo 


(2.3.31) 


because  the  error  from  each  increment  makes  the  equation  of  equilibrium  not  exactly 
satisfied,  adding  the  residual  force  to  the  next  time  step  can  not  only  reduce  the  error  but 
also  reduce  iterations  [27]  .  Finally  the  simultaneous  equations  are  in  the  form  as 
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A"J  [S"]r[D"][S"]</n  j{i«" 


=  /  [«*1  '[‘1 1  +  v“)[0"]{As;}  +  (k’ +  1){»")1  da  +  r'  (2.3.32) 


Simply 


M  (A«n}  =  {Ar} 


(2.3.33) 


where  [/if?]  is  termed  as  the  tangential  stiffness  matrix. 


[*r]  =  kn\  iBnY{DnJ_Bn-]dQ. 


(2.3.34) 


the  {A F")  is  called  pseudo  force, 


un = J  EB"]rc*"(i + + (*"  +  / 


(2.3.35) 


By  assembling  the  simultaneous  equations,  applying  the  boundary  condition  and  solving, 
we  get  {Aw*}  and  then  calculate  the  new  coordinates 


(.v"+')  =  (.0  +  {A«“) 


(2.3.36) 
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new  [5*']  is  then  modified  by  using  {.v',+I}  =  {x"}  +  (Am"}  to  form  new  stiffness  matrix  for 
iterations.  After  error  between  two  iterations  becomes  smaller  than  a  critical  value  for 
convergence,  the  same  routine  is  repeated  for  the  (n  +  1  )th  time  step. 

Solution  Technique.  There  are  several  methods  for  solving  nonlinear  simultaneous 
equations  [27]  in  general,  but  no  one  is  perfect  for  every  problem.  In  order  to  get  the 
correct  result,  the  intervals  of  time  steps  should  be  sufficiently  small  and  iterations  large 
enough  for  getting  a  convergent  and  stable  solution  at  that  step.  On  the  other  hand, 
increment  lengths  are  sought  to  be  large  and  iterations  to  be  fairly  few  to  reduce  the 
computer  time.  The  strategy  of  balancing  the  increment  intervals  and  iterations  is  a 
matter  of  experience.  A  combination  of  Euler's  rule  [27]  and  residual  force  correction 
[29]  is  used  in  this  code  to  reduce  iterations  on  each  step.  And  the  assumption  that  the 
stress  varies  linearly  during  each  time  step  is  used  so  that  the  increments  can  be  taken 
relatively  larger.  The  procedures  for  reducing  iterations  and  increasing  the  allowable 
increment  lengths  have  made  this  code  fairly  efficient  so  that  we  can  simulate  a  longer 
sintering  process 

Free  strain  rate  and  Constitutive  Parameters.  The  sintering  strain  loads  are 
derived  from  experiments  made  by  Rahaman  [1]  and  some  of  his  unpublished  data.  In 
the  process  of  isothermal  sintering,  the  temperature  is  quickly  increased  to  and  then  is 
kept  constant  at  750  °C.  The  relative  density  of  the  unreinforced  matrix  is  measured 
versus  time.  The  sintering  matrix  is  considered  to  be  viscoelastic.  Maxwell  model  was 
used  to  simulate  the  behavior  of  the  matrix  under  stresses.  The  viscosity  coefficient  ij  is 
measured  by  the  technique  of  loading  dilatometry  [15]  .  The  Young's  modulus  E  and 
Poisson's  ratio  are  assumed  [23]  . 
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Features  compared  with  NASTRAN.  This  computer  code  not  only  has  some 
viscoelastic  (creep)  analysis  functions  similar  to  that  NASTRAN  does,  but  also  has  some 
special  functions  for  the  sintering  problem  that  NASTRAN  does  not  have. 

The  following  are  the  major  differences  of  this  code  from  NASTRAN: 

In  this  code,  the  Young's  modulus  E,  Poisson's  ratio  v  and  viscosity  7  can  be  made 
to  depend  on  local  density  (not  only  overall  density)  while  in  NASTRAN  these  parametes 
can  only  be  made  to  depend  on  time  (or  temperature)  and  should  be  prescribed  apriori. 
Further,  coefficient  of  viscosity  7  can  only  be  made  stress  dependent  in  NASTRAN. 

Coefficient  of  viscosity  specifications  can  be  made  as  many  as  needed  up  to  the  total 
number  of  elements  while  in  NASTRAN  all  elements  must  have  the  same  7.  This  makes 
it  impossible  to  use  NASTRAN,  for  example,  with  the  the  three  region  composite  sphere 
model  (see  Appendix). 

This  code  can  be  run  on  supercomputers,  such  as  CNSF  at  Cornell,  New  York  and 
NSCA  at  Champaign,  Illinois,  to  simulate  the  entire  process  of  sintering.  However 
NASTRAN  as  a  commercial  package  program  is  currently  not  supported  by  some  of  these 
supercomputers.  Nevertheless,  this  computer  code  is  a  base  into  which  new  functions  can 
be  easily  implemented. 

2.3.5  Preliminary  Results  and  Discussion. 

This  finite  element  program  has  just  been  developed  and  is  still  being  tested  and 
undergoing  improvements,  only  a  few  of  simple  cases  have  been  studied  and  the  results 
for  these  cases  are  presented. 
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Single  Inclusion  Case.  A  simple  mesh  was  used  in  single  inclusion  case  (10% 


volume  fraction)  as  shown  in  Fig.  29.  The  average  relative  density  of  the  constrained 
matrix  against  the  unconstrained  matrix  was  calculated  and  compared  with  the 
experimental  data  as  shown  in  Fig.  30.  The  result  reconfirms  Scherer's  theory  as 
mentioned  earlier. 

Multiple  Inclusion  Case.  The  multiple  inclusion  model  seems  more  realistic  than 
single  inclusion  model  inspite  of  idealization  of  uniform  size  and  distribution  of  inclusions 
in  the  matrix.  Due  to  symmetry,  only  the  shaded  region  is  considered  as  shown  in  Fig. 
31.  A  rough  element  mesh  is  used  as  shown  in  Fig.  32.  A  rigid  element  (its  Young's 
modulus  is  relatively  much  larger  than  the  rest)  is  located  away  from  the  first  inclusion 
at  the  origin.  The  size  of  the  inclusions  and  distance  between  them  are  arranged  to  match 
the  inclusion  volume  fraction  and  8-noded  isoparametric  elements  are  used.  The 
difference  in  the  predicted  relative  densities  for  the  two  cases,  single  and  multiple 
inclusions,  is  very  minimal,  as  shown  in  Fig.  33.  This  is  probably  because  the  region 
considered  was  too  small  and  mesh  too  coarse. 

A  interesting  phenomenon  from  <he  preliminary  result  is  that  the  variation  of 
density  is  not  sensitive  to  the  viscosity  coefficient  but  quiet  so  to  the  poisson's  ratio.  It 
agrees  with  results  from  the  viscoelastic  analysis  (see  Appendix)  that  no  volumetric  strain 
except  elastic  part  exsits  if  the  poisson's  ratio  is  constant  throughout  the  process. 

Future  Work.  Since  this  computer  code  of  finite  element  analysis  has  just  been  set 
up,  it  needs  to  be  tested  further  against  available  closed  form  solutions  and  phvsical 
interpretation  of  the  results  to  ensure  the  validity  of  the  analyses.  This  computer  code 
will  then  be  used  in  analysing  and  investigating  the  following  problems:  The  interactions 
among  inclusions;  the  crack-like  damages  including  pre-existing  cracks  from  cold 
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compaction  and  cracks  induced  while  sintering  and  the  effect  of  matrix  damages  on  the 
densification  rates.  In  addition,  a  major  effort  will  be  launched  to  combine  non-mechanical 
factors,  such  as  mass  diffusions,  grain  surface  energy  minimization  and  grain  growth. 
Although  most  of  the  work  will  be  done  by  two  dimensional  simulations,  some  special 
problems  will  be  analysed  in  three  dimension  to  estimate  the  errors  produced  by  two 
dimensional  approximation. 
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Fig.  32.  A  Simple  Finite  Element  Mesh  of  Multiple  Inclusion  Model 


Unconstrained  sintering  (experiment) 
Constrained  with  10%  inclusions  (experiment) 


oo  n  <o  m 

d  o  o  d 


fU;su3Q  anuoidtj 


66 


Fig.  33  Predictions  for  matrix  with  10%  volume  fraction  inclusions. 


2.4  SUMMARY  OF  ACCOMPLISHMENTS 


Experimentally,  the  effect  of  key  material  and  processing  parameters  (e.g.  inclusion 
volume  fraction,  inclusion  size,  inclusion  distribution,  matrix  uniformity  and  sintering 
temperature)  on  the  sintering  of  a  model  polycrystalline  ceramic  composite  has  been 
determined.  The  composite  consisted  of  a  fine-grained  ZnO  matrix  and  an  inert,  rigid, 
particulate  inclusion  phase  of  Zr02- 

The  sintering  behavior  of  the  composite  was  shown  to  depend  significantly  on  the 
method  used  to  form  the  green  composite  body. 

It  was  demonstrated  that  ceramic  particulate  composites  with  an  appreciable  inclusion 
content  can  be  sintered  conventionally  to  nearly  full  density  if  forming  methods  which 
minimize  inclusion  particle  interactions  and  maximize  the  matrix  powder  packing 
immediately  surrounding  the  inclusions  were  used.  As  an  example,  composites  containing  up 
to  35  v%  inclusions  that  were  formed  by  coating  the  inclusion  particles  with  the  matrix 
phase  by  chemical  precipitation  followed  by  die-pressing  were  sintered  to  nearly  full 
density  under  identical  sintering  conditions  used  for  the  unreinforced  matrix. 

Improved  sintering  conditions  for  ceramic  composites  were  shown  t''  prevail  at  higher 
temperatures. 

For  composites  formed  by  conventional  processing  (i.e.  conventional  mixing  and 
die-pressing),  the  densification  rate  decreased  drastically  with  increasing  inclusion 
volume  fraction  or  with  decreasing  inclusion  size  and  increased  with  increasing 
temperature. 

Composites  formed  by  slip  casting  showed  somewhat  similar  sintering  behavior  to  that 
for  conventionally  processed  composites. 

The  main  factors  controlling  the  sintering  of  conventionally  processed  composites  and 
slip  cast  composites  were  shown  to  involve  the  packing  of  the  matrix  immediately 
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surrounding  the  inclusions  and  inclusion  particle  interactions  that  constrained  the 
matrix. 

The  effect  of  the  inclusion  phase  on  the  composite  viscosity  was  determined  for  a 
system  consisting  of  a  soda-lime  glass  matrix  containing  spherical  nickel  inclusions.  For 
10  v%  inclusions,  the  composite  viscosity  was  approximately  the  same  as  the  unreinforced 
glass  and  there  was  almost  no  effect  of  the  inclusion  size.  For  20  v%  inclusions,  the 
composite  viscosity  increased  significantly  with  decrease  in  the  inclusion  size.  The 
results  indicate  that  interaction  between  the  inclusion  particles  is  an  important  factor 
for  glass  matrix  composites  containing  more  than  *  10  v%  inclusions. 

Theoretically,  modelling  studies  were  initiated  to  account  for  the  experimental  data. 
A  finite  element  program  was  developed  to  describe  the  data  over  the  entire  sintering 
process.  The  sintering  problem  was  considered  as  a  special  thermomechanical  problem  in 
which  the  instantaneous  free  strain  rate  as  well  as  the  material  properties  are  density 
dependent.  The  effects  of  single  inclusion  and  multiple  inclusion  cases  and  inhomogeneity 
in  the  matrix  density  was  determined. 
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2-5  FUTURE  WORK 


The  experiments  performed  in  the  first  year  yielded  a  number  of  important  results.  It 
is  important  to  determine  whether  the  results  are  specific  to  the  Zn0/Zr02  system  studied. 
The  experiment  ^  therefore  be  extended  to  a  composite  system  consisting  of  a 
fine-grained  AI2O3  matrix  reinforced  with  Z1O2  inclusions. 

The  experiments  will  be  extended  to  include  whisker  reinforced  composites.  A  system 
consisting  of  a  fine-grained  AI2O3  matrix  reinforced  with  SiC  whiskers  will  be  utilized. 

It  is  important  to  determine  whether  the  effect  of  an  inclusion  phase  on  the  viscosity 
of  a  glass  matrix  is  qualitatively  different  from  that  for  a  polyciystalline  matrix.  The 
first-year  experiments  for  glass  matrix  composites  will  therfore  be  extended  to  a 
polycrystalline  matrix  (ZnO)  containing  ZrC>2  inclusions. 

Theoretically,  the  computer  code  for  the  finite  element  analysis  that  was  developed  in 
the  first  year  will  be  tested  further  against  available  closed  form  solutions  and  the 
validity  of  the  results  will  be  checked.  The  computer  code  will  then  be  used  to  analyze 
the  effects  of  multiple  inclusions,  inclusion  interactions,  and  crack-like  damage  arising 
from  compaction  and  sintering  on  the  composite  densification  rate.  The  incorporation  of 
diffusion-controlled  mechanisms  for  densification  and  grain  growth  into  the  model  will 
also  be  investigated. 


69 


.U  l 


2.6  REFERENCES 


1.  L.  C.  De  Jonghe,  M.  N.  Rahaman,  and  C.  H.  Hsuech,  'Transient  stresses  in 

Bimodal  Compacts  During  sintering",  Acta  metall.,  34  (7)  1467-1471  (1986). 

2.  R.  K.  Bordia  and  R.  Raj,  "Sintering  of  TiOz  —  AlzO}  Composites:  A  model 

Experimental  Investigation",  J.  Am.  Ceram.  Soc.,  71  (4)  302-310  (1988). 

3.  VV.  H.  Tuan,  E.  Gilbart,  and  R.  J.  Brook,  "Sintering  of  Heterogeneous  Ceramic 

Compacts,  Part  I:  A1Z03  -  AlzO}  ",  J.  Mater.  Sci.,  24  1062-1068  (1989). 

4.  M.  N.  Rahaman  and  L.  C.  De  Jonghe,  "Effect  of  Rigid  Inclusions  on  the  Sintering 
of  Glass  Ponder  Compacts,  J.  Am.  Ceram.  Soc.,  70  (12)  c-308-c351  ( 1987). 

G.  W.  Scherer,  "Sintering  with  Rigid  Inclusions",  J.  Am.  Ceram.  Soc.,  70  (10) 
719-725  (1987). 

6.  R.  Raj  and  R.  K.  Bordia,  "Sintering  of  Bimodal  Powder  Compacts",  Acta  Metall, 

32  (7)  1003-1019  (1984). 

7.  C.  H.  Hsueh,  A.  G.  Evans,  R.  M.  Cannon,  and  R.  J.  Brook,  "Viscoelastic  Stresses 

and  Sintering  Damage  in  Heterogeneous  Powder  Compacts",  Acta  Metall.,  34 
(5)  927-936  (1986). 

8.  F.  F.  Lange,  "Densification  of  Powder  Rings  Constrained  by  Dense  Cylindrical 

Cores",  Acta  Metall.,  37  (2)  697-704  (1989). 

9.  F,  F.  Lange,  "Constrained  Network  Model  for  Predicting  Densification  Behavior  of 

Composite  Powders",  J.  Mater.  Res.,  2  (1)  59-65  (1987). 


70 


10.  M.  W.  Weiser  and  L.  C.  De  Jonghe,  "Inclusion  Size  and  Sintering  of  Composite 

Powders",  J.  Am.  Ceram.  Soc.,  71  (3)  C- 125- 127  (1988). 

11.  L.  C.  De  Jonghe  and  M.  N.  Rahaman,  "Densification  of  Particulate  Ceramic 

Composites:  The  Role  of  Heterogeneities",  pp.  353-361  in  Mater.  Res.  Soc 
Svmp.  Proc..  Vo!  155:  Processing  of  Advanced  Ceramics  Edited  by  I.  A.  Aksay. 
The  Materials  Research  Society,  Pittsburg,  PA  (1989). 

12.  R.  K.  Bordia,  and  G.  W.  Scherer,  "On  Constrained  Sintering  —III.  Rigid  Inclusions" 

Acta  Metall.,  36  (10)  2411-2416  (1988). 

13.  L.  C.  De  Jonghe  and  M.  N.  Rahaman,  "Sintering  Stress  of  Homogeneous  and 

Heterogeneous  Powder  Compacts",  Acta  metall.,  36  (1)  223-229  (1988). 

14.  M.  N.  Rahaman  and  L.  C.  De  Jonghe,  "Sintering  of  Particulate  Composites  Under 

a  Uniaxial  Stress",  J.  Am.  Ceram.  Soc.,  73  (3)  602-606  (1990). 

15.  M.  N.  Rahaman  and  L.  C.  De  Jonghe,  "Sintering  of  Ceramic  Particulate 

Composites:  Effect  of  Matrix  Density",  J.  Am.  Ceram.  Soc.,  74  (2)  433-536 
(1991). 

16.  L.  C.  De  Jonghe  and  M.  N.  Rahaman,  "A  Loading  Dilatometer",  Rev.  Sei.  Instrum., 

55  (12)  2007-2010  (1984). 

17.  M.  N.  Rahaman,  L.  C.  De  Jonghe,  and  R.  J.  Brook,  "Effect  of  Shear  Stress  on 

Sintering",  J.  Am.  Ceram.  Soc.,  69  (1)  53-58  (1986). 

18.  G.  W.  Scherer  and  A.  Jagota,  "Effect  of  Inclusions  on  Viscous  Sintering",  Ceramic 

Transactions,  in  Press. 


71 


19.  R.  M.  Christensen.  Theory  of  Viscoelasticity  pp.  16-20,  2nd  Edition,  Academic,  New 


York  (1982) 

20.  A.  G.  Evans,  "Considerations  of  Inhomogeneity  Effect  in  Sintering",  J.  Am  Ceram. 

Soc.,  65  (10)  497-501  (1982) 

21.  F.  F.  Lange  and  M.  Metcalf,  "Processing-Related  Fracture  Origins:  II,  Agglomerate 

Motion  and  Cracklike  Internal  Surfaces  Caused  by  Differential  Sintering",  J. 
Am.  Ceram.  Soc.,  66  (6)  398-406  (1983). 

22.  R.  K.  Bardia  and  G.  W.  Scherer,  "On  constrained  Sintering-I.  Constituive  Model 

For  A  Sintering  Body",  Acta  Metall.,  36  (10)  2393-2397  (1988). 

23.  R.  K.  Bordia  and  G.  W.  Scherer,  "On  Constrained  Sintering- 1 1.  Comparison  of 

Constitutive  Models",  Acta  metall.,  36  (10)  2399-2409  (1988). 

24.  S.  Andareson  and  I. A.  Aksay,  "Sintering  with  Rigid  Inclusions:  Pair  Interactions", 

J.  Am.  Ceram.  Soc.,  73  (1)  54-60  (1990). 

25.  G.  W.  Scherer  and  A.  Jagota,  "2-Dimemsional  Simulation  of  Sintering  with  Rigid 

Inclusion",  in  press. 

26.  Y.  Yamada  and  K.  Iwata,  "Finite  Element  Analysis  of  Thermal- Viscoelastic 

Material",  Research  of  Productions,  24  165-170  (1°72). 

27.  O.  C.  Zienkiewicz,  Finite  Element  Method  ,  pp.  451-457,  McGraw-Hill,  New 

York, (1979). 

28.  A  Jagota  and  P.  R.  Dawson,  "Simulation  of  The  Viscous  Sintering  of  Two 

Particles",  J.  Am.  Ceram.  Soc.,  73  (1)  173-177  (1990) 


72 


29.  D.  R.  J.Owen  and  E.  Hinton,  Finite  Element  in  Plasticity  ,  pp.  275-279,  Pineridge 

Press  Limited,  Swansea  U.K.,(1980) 

30.  G.  Berg  and  P.  Grau,  "Estimating  the  Viscosity  of  Viscoelastic  Solid  Based  on  a 

three  parameter  Model",  J.  Am.  Ceram.  Soc.,  72  (11)  2187-2188  (1989) 

31.  R.  M.  Christensen,  Mechanics  of  Composite  Materials  ,  pp.  31-70  John  Wiley  and 

Sons,  New  York,  (1979) 

32.  W.  N.  Findley,  J  S.  Lai  and  K.  Onaran,  Creep  and  Relaxation  of  Nolinear 

Viscoelactic  Materials  ,  pp.  64-68,  North-Holland  Publishing  Co.  (1976) 


73 


3.  PUBLICATIONS 


1.  M.  N.  Rahaman  and  L.  C.  de  Jonghe,  "Sintering  of  Ceramic  Particulate  Composites: 
Effect  of  Matrix  Density",  J.  Am.  Ceram.  Soc.,  74  [2]  433-436  (1991). 

2.  D-Y.  Jeng  and  M.  N.  Rahaman,  "Sintering  of  Amorphous  Mullite  with  Crystalline 
Inclusions",  Ceram.  Eng.  Sri.  Proc.,  in  press. 

3.  M.  N.  Rahaman  and  C-L.  Fan,  "Sintering  of  Model  Ceramic  Particulate  Comoosites,  I: 
Conventionally  Processed  Composites",  J.  Am.  Ceram.  Soc.,  in  preparation. 

4.  C-L.  Fan  and  M.  N.  Rahaman,  "Sintering  of  Model  Ceramic  Particulate  Composites, 
II:  Slip  Cast  Composites",  J.  Am.  Ceram.  Soc.,  in  preparation. 

5.  C-L.  Hu  and  M.  N.  Rahaman,  "Sintering  of  Model  Ceramic  Particulate  Composites, 
III:  Coated  Particles",  J.  Am.  Ceram.  Soc.,  in  preparation. 

6.  R.  E.  Dutton  and  M.  N.  Rahaman,  "Effect  of  Inclusions  on  the  Viscosity  of  Glass". 

J.  Mater.  Res.,  in  preparation. 

7.  H.  Wei  and  L.  R.  Dharani,  "Effect  of  Microcracking  on  the  Densification  Rate  of  a 
Ceramic  Matrix  Composite",  Intemat.  Conf.  on  Composite  Materials,  July  15-19, 1991. 
Honolulu,  Hawaii. 

8.  H.  Wei  and  L.  R.  Dharani,  "Analysis  of  Sintering  Cracks  in  a  Ceramic  with  a  Rigid 
Inclusion",  21st  Midwestern  Mechanics  Conference,  August  13-16, 1991.  Houghton,  MI. 


74 


4.  PROFESSIONAL  PERSONNEL 


Investigators: 

M.  N.  Rahaman,  Principal  Investigator. 
L.  R.  Dharani,  Co-principal  Investigator. 
R.  E.  Moore,  Senior  Investigator. 

Graduate  Students: 

R.  E.  Dutton,  Ph.D.  candidate. 

C-L.  Hu,  Ph.D.  candidate. 

H.  Wei,  Ph.D.  candidate. 

C-L.  Fan,  M.S.  candidate. 


5.  INTERACTIONS 


5.1.  Presentations 


M.  N.  Rahaman  and  D-Y.  Jeng,  "Sintering  of  Amorphous  Mullite  with  Crystalline 
Inclusions",  Am.  Ceram.  Soc.  Engineering  Division  Meeting,  January  13-16, 1991,  Cocoa 
Beach,  FL. 

M.  N.  Rahaman,  "Constrained  Sintering  of  Ceramic  Composites  and  Thin  Films",  Invited 
Seminar,  February  21, 1991,  University  of  Washington,  Seattle. 


75 


5.2.  Consultative  Function 


Consultative  visit  by  M.  N.  Rahaman,  Principal  Investigator,  to  the  University  of 
California,  Berkeley,  Lawrence  Berkeley  Laboratory,  Center  for  Advanced  Materials, 
Structural  Ceramics  Program,  to  review  areas  of  mutual  interest  in  "Sintering  of  Ceramics 
and  Composites".  February  20, 1991.  Professor  L.  C.  De  Jonghe,  Professor  of  Ceramics  and 
Staff  Senior  Scientist. 


APPENDIX 


AN  ANALYTICAL  VISCOELASTIC  SOLUTION  OF  3-REGION  COMPOSITE 

SPHERE  MODEL 

For  the  2-region  composite  sphere  in  which  the  core  is  rigid,  the  elastic  solution  is 
as  follows:  In  the  spherical  coordinate(  r,  8,  <j>),  the  displacements  is  given  by 


“r  =  «(>•) 

«<?  =  =  0 

Equation  of  equilibrium  in  terms  of  displacements  is 


d  u  ,  2  du  _  2 u_  _  n 
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for  which  the  solution  is  given  by 
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Using  the  kinematic  relations 


sr  = 


du 

dr 


t9~s4>-  r 


the  stresses  can  be  written  as 
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The  appropriate  boundary  conditions  are 


u(a)  =  0,  o{b)  =  Pb 


where  Pb  =  0  if  no  external  traction  is  applied.  The  constants  C  and  D  can  be  determined 
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where  v,  =  ( —  )3  is  the  volume  fraction  of  the  core.  So  that 
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The  mean  stress  is 
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By  the  correspondence  principle,  tlie  viscoelastic  solution  can  be  obtained. 
Denoting 


As)  =  L{/(0}  = 
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The  Laplace  transform  of  equilibrium  equation  is 
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The  general  solution  form  of  the  above  equation  is 
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Assuming  v(r)  =  constant,  then  the  Laplace  transform  of  displacement  and  stresses  are 
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The  result  shows  that  no  matter  v  .at  viscoelastic  constitutive  models  are  used,  there  is 
no  difference  in  volumetric  strain  between  the  elastic  and  viscoelastic  cases  for  the 
composite  sphere  model  under  the  condition  of  no  external  traction  and  constant  poisson's 
ratio.  However,  the  region  outside  the  composite  sphere  consists  of  a  viscoelastic  matrix 
and  dispersed  rigid  inclusions  at  the  same  volume  fraction  as  the  core.  The  outside  region 
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can  be  treated  as  a  homogeneous  medium  with  different  constitutive  parameters  and 
shrinking  rate  from  the  intermediate  region.  Therefore  it  is  reasonable  to  assume 
Pb  0.  There  are  two  ways  to  determine  Pb,  one  of  them  is  to  use  the  so  called 
self-consistent  model,  the  other  is  a  three  region  composite  sphere  model.  It  is  necessary 
to  know  the  constitutive  parameter  and  free  strain  rate  for  the  outside  region. 

Eshelby  formula  and  Hashin-Shtrikman  equations  for  shear  and  bulk  moduli  of  a 
composite  containing  dispersed  inclusions  [31]  are 
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Kc=K+Vi(K±ri 

Using  the  three  region  composite  sphere  model,  we  can  get  the  Pb  for  the  elastic  case,  that 
is 
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The  corresponding  viscoelastic  solution  for  Pb{t)  ^  0  may  be  solved  by  the 
correspondence  principle  and  the  Laplace  transform  once  the  shear  and  bulk  relaxation 
functions  or  creep  functions  are  known. 
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The  relaxation  and  creep  function  may  be  so  complicated,  however,  as  to  be 
nontransform  type.  There  are  also  several  other  limitations  in  employing  this  solution  in 
the  analysis  of  sintering  problems.  For  example,  the  poisson's  ratio  is  not  constant  and 
the  deformation  is  not  infinitesimal  through  the  sintering  process.  However  the  nature 
of  the  solution  gives  out  some  insight  into  the  problem.  For  example,  the  viscoelastic 
volumetric  strain  is  the  same  as  the  elastic  volumetric  strain  if  the  poisson's  ratio  is 
assumed  to  constant ,  even  if  the  external  traction  Pt(f)  0. 


